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Abstract. We summarize some recent results of the authors and their collaborators, regarding 
the derivation of thin elastic shell models (for shells with mid-surface of arbitrary geometry) from 
the variational theory of 3d nonlinear elasticity. We also formulate a conjecture on the form and 
validity of infinitely many limiting 2d models, each corresponding to its proper scaling range of 
the body forces in terms of the shell thickness. 
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1. Introduction 

Elastic materials exhibit qualitatively different responses to different kinematic boundary con- 
ditions or body forces. A sheet of paper may crumple under compressive forces, but it shows a 
more rigid behavior in a milder regime. A cylinder buckles in presence of axial loads. A clamped 
convex shell enjoys great resistance to bending and stretching, but if a hole is pierced into it, 
the whole structure might easily collapse. Growing tissues, such as leaves, attain non-flat elastic 
equilibrium configurations with non-zero stress, even in the absence of any external forces. 

Such observations gave rise to many interesting questions in the mathematical theory of elas- 
ticity. Its main goal is to explain these apparently different phenomena based on some common 
mathematical ground. Among others, the variational approach to the nonlinear theory has been 
very effective in rigorously deriving models pertaining to different scaling regimes of the body 
forces [9]. The strength of this approach lies in its ability to predict the appropriate model to- 
gether with the response of the plate without any a priori assumptions other than the general 
principles of 3d nonlinear elasticity. 

The purpose of this paper is to introduce some new results and conjectures on the variational 
derivation of shell theories. They can be considered as generalizations of the results in [9], justifying 
a hierarchy of theories for nonlinearly elastic plates. This hierarchy corresponds to the scaling of 

Date: July 9, 2009. 

1991 Mathematics Subject Classification. 74K20, 74B20. 

Key words and phrases, shell theories, nonlinear elasticity, Gamma convergence, calculus of variations. 

1 



2 



MARTA LEWICKA AND REZA PAKZAD 



the elastic energy in terms of thickness h, in the limit as h — ► 0. Some of the derived models were 
absent from the physics and engineering literature before. 

Acknowledgments. M.L. was partially supported by the NSF grant DMS-0707275. and by 
the Center for Nonlinear Analysis (CNA) under the NSF grants 0405343 and 0635983. R.P. was 
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2. The set-up and a glance at previously known results 

2.1. Three dimensional nonlinear elasticity and the limiting lower dimensional the- 
ories. The equations for the balance of linear momentum for the deformation u = u(t, x) G R 3 
of the reference configuration Q C R 3 of an elastic body with constant temperature and density 
read [l]: 

(2.1) d tt u - div DW(Vu) = f, 

where DW is the Piola-Kirchhoff stress tensor, / is the external body force, and the elastic energy 
density W is assumed to satisfy the following fundamental properties of frame indifference (with 
respect to the group of proper rotations 50(3)), normalization and non-degeneracy: 

VF e R 3x3 Vi? € SO (3) W(RF) = W(F), W{R) = 0, 
UJ) W(F) >c-dist 2 (F,SO(3)), 

with a uniform constant c > 0. 

The steady state solutions to (|2.ip satisfy the equilibrium equations: — d\v DWiVu) = f which, 
expressed in their weak form, yield the formal Euler-Lagrange equations for the critical points of 
the total energy functional: 

(2.3) J{u) = [ W(Vu) - [ fu, 

Jn Jn 

defined for deformation u : Q — ► R 3 . We will refer to the term E{u) = I WiVu) as the elastic 

Jn 

energy of the deformation u. 

As a first step towards understanding the dynamical problem (|2.ip it is natural to study the 
minimizers of (|2.3p . in an appropriate function space. The questions regarding existence and 
regularity of these minimizers are vastly considered in the literature. However, due to the loss 
of convexity of W, caused by the frame indifference assumption, these problems cannot be dealt 
with the usual techniques in the calculus of variations; see [1] for a review of results and open 
problems. 

One advantageous direction of research has been to restrict the attention to domains 17 which 
are thin in one or two directions, and hence practically reduce the theory to a 2d or Id problem. 
Indeed, the derivation of lower dimensional models for thin structures (such as membranes, shells, 
or beams) has been one of the fundamental questions since the beginning of research in elasticity 
|20j. The classical approach is to propose a formal asymptotic expansion for the solutions (in 
other words an Ansatz) and derive the corresponding limiting theory by considering the first 
terms of the 3d equations under this expansion [2]. The more rigorous variational approach of 
r -convergence was more recently applied by LeDret and Raoult [13] in this context, and then 
significantly furthered by Friesecke, James and Miiller [9], leading to the derivation of a hierarchy 
of limiting plate theories. Among other features, it provided a rigorous justification of convergence 
of minimizers of (12. 3D to minimizers of suitable lower dimensional limit energies. 
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2.2. T-convergence. Recall [6\ that a sequence of functionals F n : X — > [—00, +00] defined on 
a metric space X, T-converges to the limit functional F : X — > [—00, +00] whenever: 

(i) (the T-liminf inequality) For any sequence x n — > x in X, one has F(x) < lim inf n ->oo F n (x n ). 

(ii) (the T-limsup inequality) For any x € X, there exits a sequence x n (called a recovery 
sequence) converging to x € X, such that limsup^^oo F n (x n ) < F{x). 

It is straightforward that limj^oo F n (x n ) = F(x) for any recovery sequence x n — > x. 

When X is only a topological space, the definition of T-convergence involves, naturally, systems 
of neighborhoods rather than sequences. However, when the functionals F n are equi-coercive and 
X is a reflexive Banach space equipped with weak topology, one can still use (i) and (ii) above 
(for weakly converging sequences), as an equivalent version of this definition. 

A fundamental consequence of the above definition is the following. If x n is a sequence of 
approximate minimizers of F n in X: 

lim \ F n (x n ) - inf F n \ = 0, 

n— >oo X J 

and if x n — > x, then x is a minimizer of F. In turn, any recovery sequence associated to a minimizer 
of F is an approximate minimizing sequence for F n . The convergence of (a subsequence of) x n is 
usually independently established through a compactness argument. 

2.3. A glance at previously known results. Let S be a 2d surface embedded in R 3 , which is 
compact, connected, oriented, of class C ' and whose boundary dS is the union of finitely many 
(possibly none) Lipschitz curves. By n we denote the unit normal vector to S, and it : S — > S 
is the usual orthogonal projection of the tubular neighborhood onto S. 

Consider a family {S h }h>o of thin shells of thickness h around S: 

S h = { z = x + tn(x); x E S, -h/2 < t < h/2}, < h< h , 
The elastic energy per unit thickness of a deformation u € W l,2 (S h , M 3 ) is given by: 

(2.4) E\u) = ~ / WT(Vu), 

On above the properties where the stored-energy density W : M 3x3 — ► [0, 00] is assumed to satisfy 
(|2.2p and to be C 2 regular in some open neighborhood of SO (3). 

In presence of applied forces f h G L 2 (5 fc ,R 3 ), the (scaled) total energy reads: 

(2.5) J h {u) = E\u) - i / /V 

It can be shown that if the forces / scale like h a , then the elastic energy E h (u h ) at (approximate) 
minimizers u h of J h scale like /i^ 3 , where /? = aif0<a<2 and /3 = 2a — 2ifa>2. The main 
part of the analysis consists therefore of identifying the r-limit la of the energies E h (-, S h ) 
as h — > 0, for a given scaling /3 > 0. No a priori assumptions are made on the form of the 
deformations u h in this context. 

In the case when S is a subset of M? (i.e. a plate), such T-convergence was first established for 
= [12], and later [HI [9] for all /3 > 2. This last scaling regime corresponds to a rigid behavior 
of the elastic material, since the limiting admissible deformations are isometric immersions (if 
= 2) or infinitesimal isometries (if > 2) of the mid-plate S. One particular case is j3 = 4, 
where the derived limiting theory turns out to be the von Karman theory [TT]. A totally clamped 
plate exhibits a very rigid behavior already for > [5j. In case < < 5/3, the T-convergence 
was recently obtained in [4], while the regime 5/3 < < 2 remains open and is conjectured to be 
relevant for crumpling of elastic sheets |27j . 
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Much less is known in the general case when S is a surface of arbitrary geometry. The first 
result in [13] relates to scaling (3 = and models membrane shells: the limit Jo depends only on 
the stretching and shearing produced by the deformation on S. Another study [7] analyzed the 
case /3 = 2, corresponding to a flexural shell model [2j, or a geometrically nonlinear purely bending 
theory, where the only admissible deformations are isometric immersions, that is those preserving 
the metric on S (see section 2). The energy I2 depends then on the change of curvature produced 
by the deformation. 

All the above mentioned theories should be put in contrast with a large body of literature, 
devoted to derivations starting from 3d linear elasticity (see [2] and references therein). In the 
present setting one allows for large deformations, i.e. not necessarily close to a rigid motion. The 
basic assumption of the linear elasticity is not taken for granted in our context. 

3. The Kirchhoff theory for shells: = 2 and arbitrary S 
The limiting theory for (5 = 2 is precisely described in the following result: 

Theorem 3.1. [7J (a) Compactness and the T-liminf inequality. Letu h 6 W 1,2 (S h , M 3 ) be a 
sequence of deformations such that E h (u h )/h 2 is uniformly bounded. Then there exists a sequence 
c h € M. 3 such that the rescaled deformations: 

y h {x + tn) =u h (x + th/h n) - c h : S h ° — ► M 3 , 

converge (up to a subsequence) in W 1 ' 2 to y o n, where y € W 2 ' 2 (S, M 3 ) and it satisfies: 

(3.1) (Vy) T Vy = Id a.e. in S. 
Moreover: 

h{y) = [ Qa^n^-n) <liminf-^(A 
Js h ~*° h 

(b) The recovery sequence and the T-limsup inequality. Given any isometric immersion 

y € W 2 ' 2 (S, M. 3 ) satisfying 13.1}) . there exists a sequence u h <G W 1 ' 2 (S h , M. 3 ) such that the rescaled 

deformations y h (x + tn) = u h (x + th/h^n) converge to y o n in W 1 ' 2 and: 

I 2 {y) > lim sup E h {u h ). 

In the definition of the limit functional I2, the quadratic forms Q,2(%, ■) are defined as follows: 

(3.2) Q 2 (x, F tan ) = min{Q 3 (^); (F - F) tan = 0}, Q 3 (F) = D 2 W(ld)(F, F). 

The form Q3 is defined for all F £ M 3x3 , while Q2(x, •), for a given x £ S is defined on tangential 
minors F tan of such matrices. Recall that the tangent space to SO (3) at Id is so(3). As a 
consequence, both forms depend only on the symmetric parts of their arguments and are positive 
definite on the space of symmetric matrices [8i. 

The functional hiu) measures the total change of curvature (bending) induced by the defor- 
mation y of the mid-surface S. In the form of the integrand II denotes the shape operator on S, 
while II(y) is the pull back of the shape operator of the surface y(S) under y. For any orthonormal 
tangent frame T,rj G T X S there holds: 

r\ ■ Ht = 77 ■ d T n and rj ■ H(y)r = r] ■ d r N, 

where iV : S — > M 3 is the unit normal to y(S): N(x) = d T y x cLy. 
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4. The von-Karman theory for shells: j3 = 4 and arbitrary S 

For the range of scalings /3 > 2 a rigidity argument [HJ [91 [15] shows that the admissible de- 
formations u are only those which are close to a rigid motion R and whose first order term in 
the expansion of u — R with respect to h is given by RV. The displacement field V is an ele- 
ment of the class Vi of infinitesimal isometrics on S |26j . The space Vi consists of vector fields 
V G W 2 ' 2 (S,R 3 ) for whom there exists a matrix field A G W 1,2 (S, M 3x3 ) so that: 

(4.1) d T y(x) = A(z)t and A(x) T = -A{x) a.e. x G 5 Vr G T^S. 

In other terms, V is a (first order) infinitesimal isometry if the change of metric induced by the 
deformation id + eV is at most of order e 2 . 




Figure 4.1. The mid-surface 5 and its deformation. 

For /3 = 4 the T-limit turns out to be the generalization of the von Karman functional [9] to 
shells, and it consists of two terms: 

h(V, B tan ) = \f 0.2 (x, B tan - ^{A 2 ) ta ^j +^lj s Q2 C V ( A ") " An W) • 

The quadratic form Q2 is defined as in (|3.2p and A is as in (|4.ip . The second term above measures 
bending, that is the first order change in the second fundamental form of S, produced by V. The 
first term measures stretching, that is the second order change in the metric of S. It involves a 
symmetric matrix field Btan belonging to the finite strain space: 

B = cl i2 {symVw; w G W 1,2 (S,R 3 )}. 

The space B emerges as well in the context of linear elasticity and ill-inhibited surfaces |25l [TO] . 

Theorem 4.1. [15] (a) Letu h G W 1,2 (S h ,M?) be a sequence of deformations whose scaled energies 
E h (u h )/h A are uniformly bounded. Then there exist a sequence Q h G 50(3) and c h G M 3 such 
that for the normalized rescaled deformations: 

(4.2) y h (x + tn) = Q h u h (x + h/h tn) - c h : S h ° — ► M 3 

the following holds. 

(i) y h converge in W 1,2 (S h °) to tt. 

(ii) The scaled average displacements: 

I rho/2 

V h ( x ) = -f y h (x + tn)-xdt 
h J-ho/2 

converge (up to a subsequence) in W ,2 {S) to some V G Vi. 



6 



MARTA LEWICKA AND REZA PAKZAD 



(iii) The scaled strains \ symSJV h converge weakly in I? to a symmetric matrix field B tan € B. 

(iv) h(V,B tan ) <limmi h ^ l/h 4 E h (u h ). 

(b) For every V € Vi and Btan £ B there exists a sequence of deformations u h G W 1,2 (S , M 3 ) 
such that: 

(i) the rescaled deformations y h (x + tn) = u h (x + th/h^n) converge in W 1,2 (S ho ) to 7r. 

(ii) the scaled average displacements V h given above converge in W 1)2 {S) to V . 

(iii) the scaled linearized strains \ symVy' 1 converge weakly in L 2 to Btan- 

(iv) h(V,Bt an )>hmsnp h ^ l/h 4 E h (u h ). 

The special case of this theorem for plates, that is when ScK 2 , was already proved in [9]. It 
can be shown that for a flat surface, the infinitesimal isometries coincide essentially with the out- 
of-plane displacements. Also, the space B becomes then the set of all linearized strains associated 
with the in-plane displacements in S, so the functional J4 can be written directly in terms of an 
out-of-plane and an in-plane displacement. The Euler-Lagrange equations derived from this limit 
functional lead to the von-Karman equations [11] . 

5. The linear theory for shells: (3 > 4 and arbitrary S; (3 = 4 and approximately 

ROBUST S 

It was shown in [15] that for a certain class of surfaces, referred to as approximately robust 
surfaces, the limiting theory for (3 = 4 reduces to the purely linear bending functional: 

(5.1) I lin (V) = ^^Q 2 (x,(V(An)-AU) tan ) dx W G Vi, 

This class of surfaces is given by the property that any first order infinitesimal isometry V £ Vi 
can be modified to be arbitrarily close to a W 1,2 second order isometry: 

V 2 = {VEVr, (A 2 )tan£B} = V 1 . 

Convex surfaces, surfaces of revolution and developable surfaces belong to this class [15] . 

In [15], it was also proved that the T-limit of E h /h@ for the scaling regime (3 > 4 is also 
given by the functional (|5.ip . This corresponds to the linear pure bending theory derived in 
[2] from linearized elasticity. The important qualitative difference between this theory and the 
limiting theory for (3 = 4 and an approximately robust surface is in the type of convergences one 
establishes for a sequence u h satisfying E h {u h ) < Ch^ ' . Indeed, if /3 > 4, the best one can prove 
is the convergence in W 1,2 , up to a subsequence, of the rescaled displacement fields: 

I fho/2 

V h (x) = j ^ ^ y h (x + tn)-x dt 

to an element V €LV\. Note the finer rescaling parameter /i^/ 2_1 with respect to (|4.1|) . 

6. Intermediate theories for plates and convex shells: (3 e (2,4) 

In paper [FTj we focused on the range of scalings 2 < (3 < 4, looking hence for an intermediate 
theory between those corresponding to (3 = 2 and (3 > 4. On one hand, modulo a rigid motion, 
the deformation of the mid-surface must look like id + eV, up to its first order of expansion. 
On the other hand, the closer [3 is to 2, the closer the mid-surface deformation must be to 
an exact isometry of S. To overcome this apparent disparity between first order infinitesimal 
isometries and exact isometries in this context, one is immediately drawn to consider higher order 
infinitesimal isometries which lay somewhat between these two categories. This will be the subject 
of discussion in section [7l Another angle of approach, which turns out to be useful in special cases, 
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is to study conditions under which, given V G Vi, one can construct an exact isometry of the form 
id + eV + e 2 w e , with equibounded corrections w e . This is what we refer to as a matching property. 
If S C M 2 represents a plate, the above issues have been answered in [9]. In this case: 

(i) The limit displacement V must necessarily belong to the space of second order infinitesimal 
isometries: V2 = {V £ Vi; {A 2 ) tan £ where the matrix field A is as in (14.11) , 

(ii) Any Lipschitz second order isometry V G V2 satisfies the matching property. 

Combining these two facts with the density of Lipschitz second order infinitesimal isometries in V2 
for a plate [22] , one concludes through the T-convergence arguments that the limiting plate theory 
is given by the functional (15. ip over V2. Note that, for a plate, V G V2 means that there exists an 
in-plane displacement w G W l,2 {S, R 2 ) such that the change of metric due to id + eV + e 2 w is of 
order e 3 . Also, in this case, an equivalent analytic characterization for V = (V 1 ,V 2 ,V 3 ) G V2 is 
given by (V 1 ,V 2 ) = {-uiy,ujx) + (61,62) and: detV 2 U 3 = 0. 

Towards analyzing more general surfaces S, we derived a matching property and the corre- 
sponding density of isometries, for elliptic surfaces. We say that S is elliptic if its shape operator 
II is strictly positive (or strictly negative) definite up to the boundary: 

(6.1) VxeS VreT x S ^|r| 2 < (iI(x)t) • r < C\t\ 2 . 

The novelty here is the fact that for an elliptic surface, all sufficiently smooth infinitesimal isome- 
tries satisfy the matching property: 

Theorem 6.1. [T7j Let S be elliptic as in (|6.ip . homeomorphic to a disk and let for some a > 0, 
S and dS be of class C 3,a . Given V G V\ D C 2,a {S), there exists a sequence w e : S — > M 3 , 
equibounded in C 2 ' a (S), and such that for all small e > the map u e = id + eV + h?w e is an 
(exact) isometry. 

We apply this result to construct the recovery sequence in the T-limsup inequality. Clearly, 
Theorem 16. II is not sufficient for this purpose as the elements of Vi are only W 2,2 regular. In most 
T-convergence results, a key step is to prove density of suitably regular mappings in the space 
of mappings admissible for the limit problem. Results in this direction, for Sobolev spaces of 
isometries and infinitesimal isometries, have been shown and applied in the context of derivation 
of plate theories. The interested reader can refer to [24\ [2"2] for statements of these density 
theorems and their applications in [9l [3] . 

In general, even though V\ is a linear space, and assuming S to be C°°, the usual mollification 
techniques do not guarantee that elements of Vi can be approximated by smooth infinitesimal 
isometries. An interesting example, discovered by Cohn-Vossen [26], is a closed smooth surface of 
non-negative curvature for which C°° H V\ consists only of trivial fields V : S — > M 3 with constant 
gradient, whereas C 2 n V\ contains non-trivial infinitesimal isometries. Therefore C°° n V\ is not 
dense in Vi for this surface. We however have: 

Theorem 6.2. [17] Assume that S is elliptic, homeomorphic to a disk, of class C m+2 ' a up to the 
boundary and that dS is C m+ ' a , for some a G (0, 1) and an integer m > 0. Then, for every 
V G Vi there exists a sequence V n G Vi D C m,a (S,M. 3 ) such that: 

lim \\V n - V\\ W 2,2 (S ) = 0. 

Ultimately, and as a consequence of Theorems 16.11 and 16.21 the main result of [T7] states that 
for elliptic surfaces of sufficient regularity, the T-limit of the nonlinear elastic energy (12. 4j) for the 
scaling regime 2 < (3 < 4 (and hence for all (3 > 2) is still given by the functional (|5.ip over the 
linear space Vi: 
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Theorem 6.3. [17J Let S be as in Theorem\6J\ and let 2 < (3 < 4. 

(a) Assume that for a sequence of deformations u h G W 1,2 (S h , M 3 ) their scaled energies E h (u h )/h /3 
are uniformly bounded. Then there exist a sequence Q h G SO (3) and c h G R 3 such that for the 
normalized rescaled deformations in ( f^.l?| ) the following holds. 

(i) y h converge in W 1,2 (S h °) to ir. 

(ii) The scaled average displacements: 

I rho/2 
V h (x) = -^j^ j ^ ^ y h (x + tn)-x dt 

converge (up to a subsequence) in W 1,2 (S) to some V G Vi- 

(iii) I lm (V) < liminf^oV^ E h {u h ). 

(b) For every V G Vi there exists a sequence of deformations u h G W 1,2 {S h , M 3 ) such that: 

(i) the rescaled deformations y h (x + tn) = u h (x + th/hon) converge in W l,2 (S h °) to it. 

(ii) the scaled average displacements V h given above converge in W 1,2 (S) to V . 

(iii) I lm (V) > Hmsup^o l/hPE h {u h ). 

One can actually prove that for 2 < f3 < 4 and surface S of arbitrary geometry, the part (a) or 
Theorem 16.31 remains valid, and moreover 

1/^/2-1 

sym W converge (up to a subsequence) in 

L 2 (S) to 1/2{A 2 ) 

tam where A is related to V by (|4.ip . The novelty with respect to the equivalent 
result for = 4 in the first part of Theorem 14.11 is the constraint V G V-2- If S is an elliptic 
surface of sufficient regularity, the set B coincides with the whole space L 2 ym (S, R 2x2 ) [15], hence 
the constraint is automatically satisfied for all V G Vi . In the general case where S is an arbitrary 
surface, a characterization of this constraint and the exact form of B may be complicated. 

The case of the scaling range 2 < f3 < 4 is still open for general shells. The following section 
is dedicated to the presentation of a conjecture on this problem, stating that other constraints, 
similar to the inclusion V G V2, should be present for values of (3 closer to 2. Heuristically, the 
closer (3 is to 2, we expect V to be an infinitesimal isometry of higher order. 

7. A CONJECTURE ON THE INFINITE HIERARCHY OF SHELL MODELS 

If the deformations u h of S h satisfy a simplified version of Kirchhoff-Love assumption: 

u h (x + tn) = u h (x) + tN h (x), 
vector N h being the unit normal to the image surface u h (S), then formal calculations show that : 

(7.1) E\u h )^ [ \5g s \ 2 + h 2 [ \5TL S \ 2 . 

Js Js 

Here by 5g$ and 511$ we denote, respectively, the change in the metric (first fundamental form) 
and in the shape operator (second fundamental form), between the surface u h (S) and the reference 
mid-surface S. For a more rigorous treatment of this observation see e.g. [1]. The two terms in 
(|7.ip correspond to the stretching and bending energies, and the factor h 2 in the bending term 
points to the fact that a shell undergoes bending more easily than stretching. For a plate, the 
latter energy is known in the literature as the Foppl-von Karman functional. 

Another useful observation is that for minimizers u h , the energy should be distributed equi- 
partedly between the stretching and bending terms. When E h {u h ) ~ h 2 , then equating the order 
of both terms in (17. lh we obtain in the limit of h — > 0: 

[\5g s \ 2 ^0 and ^E h (u h ) « f \5TL S \ 2 . 
J S J s 
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This indeed corresponds to the Kirchhoff model, as in Theorem 13.11 [7] , where the limiting energy 
I2 is given by the bending term (measuring the change in the second fundamental form) under 
the constraint of zero stretching: 5gs = 0. As we have seen, the limiting deformation u must be 
an isometry (Vit) T Vu = Id and hence preserve the metric. 

To discuss higher energy scalings, assume that: 

(7.2) E h (u h ) ^h 13 , > 2. 

Then, as mentioned before, by the rigidity estimate [8], the restrictions of u h to S have, modulo 
appropriate rigid motions, the following expansions: 



^ = id + E A 

i=i 



Thus, 511 is of the order e and after equating the order of the bending term in (|7.1H by (|7.2jl . we 
arrive at: h 2 e 2 = far, that is: 

(7.3) e = h,M 2 - 1 . 

On the other hand, the stretching term has the form: 5gs = (Vu )Vu — Id = Yli^i e *^-«> with: 

A i= Yl sym(jy Wj ) T Vw t 

j+k=i 

indicating the i-th order change of metric. Taking into account (|7.2p . this yields: e 2i f s \Ai\ 2 ~ ', 
and so in view of (|7,3p : H^H 2 ^ ~ /j/ 3 -*^- 2 ) . A first consequence is that A\ must vanish in the 
limit as h —* 0, that is the limiting deformation is a first order infinitesimal isometry. For % > 1, 
we observe that [|-Ai||?a ~ h^ % ~ l ' l ^ i '~^\ where: 

0i=2 + -r^—. 



We conclude that if (3 < /3jv, then H^H^a ~ for i < N, and if (3 = /9jv, then ||^4tv||l 2 = 0(1). The 
study of the asymptotic behavior of the energy l/h /3 E h leads us hence to the following conjecture. 

Conjecture 7.1. The limiting theory of thin shells with midsurface S, under the elastic energy 
scaling (3 > 2 as in |7.£| ) is given by the following functional Ip below, defined on the space Vn of 
N-th order infinitesimal isometrics, where: 

The space Vn is identified with the space of N -tuples (V±, . . . , Vn) of displacements V{ : S — > M 3 
(having appropriate regularity), such that the deformations of S: 

N 

u e = id + ^2 6% 

i=l 

preserve its metric up to order e N . We have: 

(i) When (3 = f3 N +i then Ip = f s Q 2 (x, S N+1 g s ) + f s Q2 (x, 6iTl s ), where S N +igs is the 
change of metric on S of the order e^"*" 1 , generated by the family of deformations u € and 
5iUs is the corresponding first order change in the second fundamental form. 

(ii) When (3 G ((3 n +i,(3n) then Ip = J s Q 2 (x, SiU s ). 
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(iii) The constraint of N-th order infinitesimal isometry Vn may be relaxed to that of Vm, 
M < N, if S has the following matching property. For every (V\, . . . Vm) € Vm there exist 
sequences of corrections V M+1 , . . . V N , uniformly bounded in e, such that: 

M N 

u e = id + J^eV<+ Yl ^ 

i=l i=M+l 

preserve the metric on S up to order e . 

This conjecture is consistent with the so far established results in [15] for N = 1 (i.e. (3 > 02 = 4) 
and arbitrary surfaces. Note that in the case of approximately robust surfaces, any element of Vi 
can be matched with an element of V2, and hence the term J s \S2gs\ 2 in the limit energy can be 
dropped. The second order infinitesimal isometry constraint V2 is established for all surfaces when 
2 < (3 < 4. In the particular case of plates, any second order isometry in a dense subset of V2, 
can be matched with an exact isometry [9]. As a consequence, the theory reduces to minimizing 
the bending energy under the second order infinitesimal isometry constraint. A similar matching 
property for elliptic surfaces, this time for elements of Vi, is given in [T7] (see Theorems 16.11 and 
16. 2p . As a consequence, for elliptic surfaces, the limiting theory for the whole range [3 > 2 reduces 
to the linear bending. 

The case 2 < (3 < 4 remains open for all other types of surfaces. The main difficulty lies in 
obtaining the appropriate convergences and the limiting nonlinear constraints: 

£ sym((VV i ) T Vy fc ) =0, l<i<N, 

j+k=i 

for the elements of Vn when (3 < (3n- The above nonlinearity, implies a rapid loss of Sobolev 
regularity of Vi as i increases. Moreover, applying methods of [17] to surfaces changing type leads 
in this context to working with mixed- type PDEs. 
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